Cavitation from bulk viscosity in neutron stars and quark stars 
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Abstract 

The bulk viscosity in quark matter is sufficiently high to reduce the effective pressure below the corresponding vapor pressure during 
density perturbations in neutron stars and strange stars. This leads to mehanical instability where the quark matter breaks apart into 
fragments comparable to cavitation scenarios discussed for ultra-relativistic heavy-ion collisions. Similar phenomena may take place in 
kaon-condensed stellar cores. Possible applications to compact star phenomenology include a new mechanism for damping oscillations 
and instabilities, triggering of phase transitions, changes in gravitational wave signatures of binary star inspiral, and astrophysical 
formation of strangelets. At a more fundamental level it points to the possible inadequacy of a hydrodynamical treatment of these 
processes in compact stars. 



PACS numbers: 97.60.Jd, 12.38.Mh, 21.65.Qr. 26.60.Kp 

It is well-known that the bulk viscosity in dense nuclear 
matter containing hyperons as well as in quark matter can 
be very high and therefore may play an important role in 
the damping of compact star oscillations and instabilities 

fl a a i a a a i i9i iioi lui 111 Hi Q [a Hi . But so 

far it has not been realized that bulk viscosity may in fact 
be so high that a hydrodynamical treatment of the oscilla- 
tions and instabilities becomes unphysical because the ef- 
fective pressure can become lower than the corresponding 
vapor pressure and cause the fluid to break up into frag- 
ments. This may introduce interesting new phenomenol- 
ogy and at least call for a re-evaluation of various processes 
in compact stars. The effect is the high chemical potential, 
low temperature analogy of the recent demonstration at low 
chemical potential and high temperature that an expanding 
quark-gluon plasma formed in ultra-relativistic heavy-ion 
collisions may be mechanically unstable and cavitate, i.e. 
break apart into droplets 113, H [HlzSl- In the following 
cavitation and mechanical instability will be used as syn- 
onymous with such a process, and the term vapor pressure 
Pvap is borrowed from ordinary fluid theory to denote the 
pressure below which the fluid can cavitate. 

The term effective pressure Pgff will be used for the 
space-space components of the energy -momentum tensor, 
which enter the relativistic hydrodynamical equations. 
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where P is the normal (thermodynamic) fluid pressure, C, is 
the bulk viscosity, and = with fluid four- velocity 
u'' is the expansion scalar, which has dimension of 1/time 
and can be thought of as the rate of expansion (see below). 
Only homogeneous and isotropic perturbations will be dis- 
cussed. Cavitation will happen when Pgff becomes smaller 
than the corresponding vapor pressure, Pyap- In the case 
of heavy-ion collisions or self-bound strange quark matter 
Pvap is the vacuum pressure, i.e. Pyap = 0, but in the case 
of a phase boundary inside a compact star Pyap = Pph, i-e. 
the pressure at the phase boundary. Terms including shear 
viscosity are neglected; such terms could further reduce 



Peff, but in most parameter ranges of interest for cavitation 
in compact stars the effects of bulk viscosity will dominate. 
Only first-order (Navier-Stokes) contributions to pressure 
are considered since it has been found that inclusion of 
second-order terms do not introduce major changes at least 
in the context of heavy-ion physics [19]. The bulk viscosity 
calculated from lattice QCD at the high temperatures and 
low chemical potentials relevant for heavy-ion collisions 
is at the border-line of being large enough to drive the ef- 
fective pressure negative on the short (strong interaction) 
time-scales involved iflTl (for 1 -dimensional 

boost expansion B = where t is the proper time). 
The same lattice calculations are relevant for the cosmo- 
logical quark-hadron transition. Here = 2>H where H 
is the Hubble-parameter. For H = l/2t and t ~ 10^^ s 
the effects of the bulk viscosity in the cosmological quark- 
hadron transition are negligible. 



But there is another astrophysical setting with a much 
higher bulk viscosity where mechanical instability could 
be relevant, namely neutron stars or quark stars. For ex- 
ample Peff < occurs for a pressure of 10'^^ dyn/cm-^, 
a bulk viscosity of 10'^^ g/cms, and an expansion scalar 
G > 10'' s^^, numbers which are realistic in a compact star 
setting. In the following the conditions necessary for cav- 
itation are derived, the relevant adiabatic indices for quark 
matter in compact stars are calculated, and the parameter 
ranges where Pgff < Pyap are presented. Mechanical insta- 
bility turns out not to occur for nuclear matter, but it can 
play an important role in compact stars containing kaon 
condensation or quark matter. This may have interesting 
consequences for compact star phenomenology. 

The conditions necessary for mechanical instability in 
dense matter are best described in terms of adiabatic in- 
dices, defined as F = d In P/ d In n^, where P is the pres- 
sure and ris is the baryon density. As demonstrated in jvl] 
the bulk viscosity in a dense, relativistic fluid can be ex- 
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where iv is the frequency of an oscillation, r is the equili- 
bration time-scale for the relevant microscopic processes, 
and and F^q are the adiabatic indices for frozen chemi- 
cal composition and equilibrium composition respectively. 
The bulk viscosity measures the impact of deviation from 
equilibrium via the competition between the time-scale 
characterizing the microscopic response, r, and the vari- 
ation in fluid density with frequency to and it enters the 
relativistic hydrodynamical equations in the combination 
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This equation demonstrates that the effect of bulk viscosity 
is largest for lot ~ 1, and it shows that the condition for 
cavitation, Peff < -fvap> is 

(Fpr - rEq)a;r G . P 
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In particular a necessary condition for cavitation is 



A = (Fpr - T 
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for which Peff < ^vap in the following range. 



A - JA^ - A2 A + JA^ - A2 



For a sinusoidal volume oscillation with amplitude A 
(V = Vo(l + Asm{ujt))) and non-relativistic flow ve- 
locities, B = Aiv cos{ujt) /{I + A sm{ujt)), so for not too 
large amplitude Ac ~ 2 (1 — Pmp/P) / Acos{ijjt). This 
shows that cavitation is most likely near a phase bound- 
ary {P —5- i-*vap)> that it is facilitated by a large ampli- 
tude, and that the effect of bulk viscosity is most important 
when cos(u;t) — > 1, i.e. when the expansion rate is maxi- 
mal, whereas the impact on phase-transformation from the 
change in thermodynamic pressure is most important for 
sin(a;t) —1. 

The bulk viscosity of nuclear matter in neutron stars with 
and without hyperons has been presented for various equa- 
tions of state in the literature Id S II, [HI [HI]- Published 
numbers are often close to the interesting regime for me- 
chanical instability, but the necessary condition Eq. ([5]) may 
be hard to fulfill unless the amplitude is very large. For 
example the values of Fpr and Tpq published in Ref. 112 ill 
differ at most by 0.8. The highest published value for A 
(for neutron stars with kaon condensates Bllil ^ is around 5. 
For equations of state that involve phase transitions, such 
as neutron stars with a kaon condensed core ll22ll . the rele- 
vant vapor pressure for the core is the pressure of the low- 
density phase at the core boundary, so Ac near the 



boundary, which means that there will be an outer layer in 
the core where cavitation triggered by bulk viscosity can 
take place for a range of ujt in addition to phase trans- 
formation caused by the oscillation in the thermodynamic 
pressure. 

Bulk viscosities in dense quark matter are known to be 
even higher than in nuclear matter for some ranges of pa- 
rameters fEllSSMIiEIillilHIil, and there- 
fore the most likely setting for cavitation is quark matter in 
strange stars if quark matter is absolutely stable or in the 
inner parts of neutron stars if quark matter becomes stable 
only at high pressure. As demonstrated below, the neces- 
sary condition for mechanical instability is indeed fulfilled 
in the outer tens of meters of strange stars, and for other 
parameters where strange matter is not absolutely stable 
cavitation can occur in quark matter cores deeper inside 
neutron stars. 

Quark matter will be treated within the MIT bag model 
where the pressure is given by 

P=-B + PM + PdiP^d) + Psil^s) + Peif-ie), (7) 

where B is the bag constant to be thought of as the en- 
ergy density of the confining vacuum, and Pi{ni) are 
the Fermi-gas contributions from the constituent particles, 
which are up, down and strange quarks as well as elec- 



trons 0231 124 125L 12611 . The star is assumed transparent to 
neutrinos. The temperature is set to zero for calculations 
of adiabatic indices (but not for calculations of r), which 
is an excellent approximation after the first few seconds 
in the life of a proto-neutron star. The pressure contribu- 
tions are (for massless up and down quarks and electrons, 
and strange quarks with mass m and Fermi momentum 

k = (,,2 _ ^2^1/2^ 
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and the corresponding number densities are ni{p,i) = 

dP/dfi,. 

The total baryon density Ub = 
[nuifJ-u] + nai^J-d) + risiUs)] and local charge 
neutrality (which is obeyed under most conditions except 
in a mixed quark-hadron phase) gives the constraint 

If reactions are very slow compared to the time-scale 
for density change, the composition is effectively frozen, 
and the highest possible value of the adiabatic index is 
obtained. To calculate Fpr = din P/d Inns Ipreeze = 
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Fixed composition means that the individual densi- 
ties vary in proportion to the baryon density, leaving 
Xi = Ui/nB constant, so for example dfiu/dus = 
d{ir'^XunBy^^/dnB = jiu/'inB, and dfis/dnB = 
k^/SfisTT-B- Collecting terms this gives 

4[P + 5 + (n,A;2/4/i,-P,)] 



Fr 



3P 
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reproducing the well-known limit for massless quarks, T = 
4(P + B)/3P, which approaches the equally well-known 
value r = 4/3 for a gas of extremely relativistic particles 
in the limit where P ^ B. 

If all particle reactions are very fast compared to the dy- 
namical time-scale for density change, the composition has 
time to equilibrate, and the adiabatic index will take its 
minimum value, r^q. The values of Xi will no longer be 
constant. Instead chemical equilibrium conditions are de- 
termined by the possible processes, which are 

s + <-> n + d (12) 

d <-> n + + Fg 



S ^ U + e + Ue 

leading to constraints on the chemical potentials 



(13) 



which together with the definition of baryon density and 
the condition for charge neutrality gives 4 constraints on 
the 5 variables taken to be the four chemical potentials plus 
Ub or the four particle densities and Ub, which again al- 
lows the calculation of all necessary derivatives in Eq. 
as a function of n^. The calculation involves implicit dif- 
ferentiation of a fairly long expression for the relation be- 
tween n„ and Ub ■ Therefore only the numerical results are 
given here. Again, in the limit of massless quarks, the ex- 
pression for converges to 4(P + B)/3P, which is a 
useful consistency check on the results, and also confirms 
that Fgq = Fpr for massless quarks. 

To appreciate the numerical results for A = (Fpr — '^Eq) 
shown in Figure 1 as a function of baryon density, Ub, it 
is important to note that the surface of a stable, self-bound 
strange star (possible for 145 MeV < B'^^'^ < 165 Me V, 
depending on the value of m lE3[l24ll25il23l ) is character- 
ized by P = 0, and that the density increases from zero to 
well above nuclear matter density when crossing the sur- 
face. For massless quarks the surface density is exactly 4S 
and contrary to an ordinary neutron star (where the density, 
not just the pressure, goes to zero at the stellar surface), 
this density increases only by a factor of a few from stel- 
lar surface to stellar center 1231 l24ll25L 12611 . This unusual 
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FIG. 1: Difference in adiabatic indices, A = Fpr — Tpq, plot- 
ted against baryon number density, ub for bag constant B^^^ = 
155 MeV and strange quark mass m as indicated. All curves are 
diverging near the density where P — (which varies slightly 
with m for fixed B). Except for a shift to higher density, similar 
curves are obtained for higher values of B where strange matter 
is unstable but may form a core in neutron stars. 



behavior of self-bound matter near the stellar surface is re- 
flected in a very special behavior of the adiabatic indices. 
Again for massless quarks one sees that F = 4{P+B)/3P 
diverges at the surface, since P ^ 0, so non-interacting 
quark matter has the special property that the equation of 
state is extremely stiff (high F) at the lowest (albeit still 
very high) densities, whereas the equation of state is soft 
(F 4/3) at the highest densities. And for m > it 
leaves room for a significant numerical difference between 
the adiabatic index calculated for frozen composition and 
equilibrium composition respectively as illustrated in Fig- 
ure 1. 

Figure 1 illustrates that the necessary condition for me- 
chanical instability, Eq. ([5]), is obeyed near the surface of a 
strange star, where cavitation involves Pyap = and there- 
fore Ac ~ 2/^cos(u;t) for sinusoidal perturbations. For 
a typical compact star mass of 1.4 solar masses A > 10 
in a surface layer stretching to a depth of 100 meters and 
A > 100 in the outer few meters. If quark matter is not 
absolutely stable but found only as a quark matter core in- 
side a neutron (or so-called hybrid) star, the outer meters 
of this core will obey Eq. dS]) for high A. In this case the 
relevant vapor pressure is the pressure of the nuclear mat- 



ter phase at the core boundary, P^ 





near the quark core boundary where P — > Pph. If quark 



vap 



Pph, so Ac 
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matter occurs in a mixed quark-hadron phase the results 
above cannot be directly applied, because quark matter in 
such a mixed phase does not obey local charge neutrality. 
However there is reason to believe that a similar effect will 
occur there. 

When the necessary condition for cavitation is fulfilled, 
the actual behavior is determined by ujt. The microscopic 
relaxation time, r, has been calculated for the relevant pro- 
cesses, and can take essentially all values from 10~^° — 10 s 
depending on temperature and strange quark mass in par- 



ticular 

llirilllllEliEIillillllllll- Cavitation is 
most likely for ujt ~ 1 (Eq. and this can very well 
occur, since many dynamical processes in compact stars 
reach such time-scales, including stellar radial oscillations, 
r-mode instabilities |27], and tidal deformation during the 
final stages of binary compact star inspiral [28]. Therefore 
the consequences can be important for detailed studies of 
these phenomena and related observables such as the grav- 
itational wave emission signatures, pulsar glitches and pos- 
sible triggering of phase transitions that could be relevant 
for gamma-ray bursts and other energetic sources. Also, 
the detailed break-up process in binary inspiral of strange 
stars will decide the mass distribution of fragments in form 
of strangelets sought for in cosmic rays [28, 29, 30]. 

Quark matter has here been treated within the simplest 
version of the MIT bag model without inclusion of gluon- 
exchange corrections. Such corrections have been shown 
(as far as the equation of state is concerned for small strong 
coupling) to be equivalent to changing the value of B [l24ll . 
so the qualitative scenario should remain unchanged. Go- 
ing beyond the bag model will clearly change the numeri- 
cal estimates, but the possibility of having large and even 
divergent values of the adiabatic indices, and therefore pos- 
sibly also large values of A, appears linked to the spe- 
cial properties of self-bound matter rather than to the spe- 
cific equation of state. Quark pairing as is expected to oc- 
cur with color superconductivity and color-flavor locking 
in the infinite density limit [31] may have only minor ef- 
fect on the equation of state and adiabatic indices, but it 
significantly (often exponentially) decreases the bulk vis- 
cosity via increasing the characteristic microscopic time- 
scale, rllSIilEIilElIil-This would significantly 
change the mechanical instability window, at least for full 
color-flavor locking whereas the change in bulk viscosity is 
smaller for phases with partial pairing, such as 2SC. How- 
ever, pairing is expected to be most important at the highest 
densities, and it could well be that quark matter is unpaired 
at the lowest densities where the mechanical instability oc- 
curs. 

The bulk viscosity in dense quark matter has been shown 
to make the effective pressure smaller than the vapor pres- 
sure under conditions relevant for neutron stars or strange 
stars. This indicates the onset of cavitation where the quark 
matter breaks apart into fragments which may survive as 
strangelets (if quark matter is stable) or hadronize in a 
manner comparable to recent scenarios for ultra-relativistic 



heavy-ion collisions. At the very least it points to the 
need for a more detailed investigation of instabilities and 
oscillations in compact stars containing quark matter. In 
particular, the cavitation process represents a new mecha- 
nism for damping oscillations, r- modes, and other pertur- 
bations in compact stars. Similar phenomena can occur in 
the outer parts of kaon condensed cores in neutron stars 
(and for other systems with a phase transition separating a 
core from the outer parts). Other types of dense nuclear 
matter do not seem to fulfill the necessary condition for 
cavitation, but the numbers are sufficiently close that also 
this may deserve further study. The effect of bulk viscosity 
is complementary to phase transformations caused by the 
perturbation in the thermodynamic pressure itself. The in- 
terplay between these contributions and questions related 
to the spectrum and growth rates of cavities lIlSll remain to 
be studied. 
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